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Abstract Collective flow, its anisotropies and its event-to-event fluctuations in relativistic 
^ . heavy-ion collisions, and the extraction of the specific shear viscosity of quark-gluon plasma 

(QGP) from collective flow data collected in heavy-ion collision experiments at RHIC and LHC 
qq ■ are reviewed. Specific emphasis is placed on the similarities between the Big Bang of our universe 

£NJ , and the Little Bangs created in heavy-ion collisions. 
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1 HISTORICAL OVERVIEW 

The idea that high energy collisions between hadrons and atomic nuclei can be 
described macroscopically using relativistic fluid dynamics has been around for 
over half a century ([TJ) . It gained strength in the late 1970's to mid 1990's as 
fixed-target heavy-ion collision experiments at beam energies from a few hundred 
MeV to 200 GeV per nucleon revealed unmistakable evidence for the formation of 
dense matter that underwent collective expansion in the directions perpendicular 
to the beam ("transverse flow") ([2|). The discovery of transverse, in particular 
elliptic flow (anisotropic emission around the beam direction) led to the the- 
oretical development of numerical codes that solved the equations of relativistic 
ideal fluid dynamics in one and two transverse directions, using coordinates and 
initial conditions appropriate for relativistic heavy- ion collisions ([5]). The pre- 
dictions from such models for the momentum distributions of the finally emitted 
hadrons reproduced qualitatively all experimentally observed features of soft (i.e. 
low transverse momentum) particle production in heavy- ion collisions ([6]), but 
failed in important aspects quantitatively (for example, hydrodynamics overpre- 
dicted the observed elliptic flow by about 50% (JTj) at CERN SPS energies while 
microscopic approaches based on a kinetic description of systems of scattering 
hadrons got it about right but seriously underpredicted it later at higher 
energies ([9])). 

In 2000 the Relativistic Heavy Ion Collider (RHIC) began operation at ten 
times larger center of mass energy than previously available. It was built to 
mass-produce quark-gluon plasma (QGP), a new type of dense matter in which 
hadrons dissolve into deconfined colored degrees of freedom, quarks and gluons, 
that was predicted by the modern theory of the strong interaction, Quantum 
Chromodynamics (QCD). RHIC data appeared to agree, for the first time, quan- 
titatively with predictions from ideal fluid dynamics (JSj) • This caused a paradigm 
change in the field: Having expected gas-like behavior, based on the ideas of 
asymptotic freedom and color Debye screening in the QGP, the community was 
forced to accept instead the notion that the QGP is a strongly-coupled plasma 
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that flows like a liquid (flOlfTT]) . 

Later it became clear that some of the initial success of ideal fluid dynamics 
was artificial and due to the use of an inadequate equation of state for the fireball 
matter (|12p as well as an incorrect treatment of the chemical composition of the 
fireball during its late hadronic stage (|13p . A key insight was (I14p i5 p6l) that the 
hydrodynamic description should be restricted to the dense quark-gluon plasma 
(QGP) stage of the collision whereas the late hadronic stage (after the recombina- 
tion of colored quarks and gluons into color-neutral hadrons) is too dissipative for 
a fluid dynamical approach (|17p and must be described microscopically. Hybrid 
codes that couple (ideal) fluid dynamics of the QGP to a microscopic hadron 
cascade worked better than a purely hydrodynamic approach (|14|, I15 ^ [T6|) and 
explained why pure fluid dynamics was quantitatively successful only in central 
collisions between large (A ~ 200) nuclei at midrapidity at top RHIC energy, but 
gradually broke down in smaller collision systems, in more peripheral collisions, 
away from midrapidity, and at lower collision energies (|18p . 

This observation focused the attention of the community on the importance of 
dissipative effects which had thus far been ignored. With improving experimental 
precision it became increasingly obvious that even in the most central collisions 
among the heaviest available nuclei the data required some degree of QGP viscos- 
ity (|16y i9ir20). The basic success of the ideal fluid approach suggested, however, 
that this viscosity had to be small. This begged the question "How small ex- 
actly?" which marked the transition of the heavy-ion program from its initial 
QGP discovery phase to a second stage of quantitative characterization of the 
QGP. 

Theoretical work on strongly coupled quantum field theories, borrowing tools 
from superstring theory, established a lower limit around for the specific 

shear viscosity n/s of the QGP (|2ip (see Chpt.[2|). How to extract this quantity 
from experimental data will be described in this review. In this extraction a key 
role is played by the fact that the initial conditions in heavy-ion collisions fluctuate 
from event to event (f22| l2~3, 24. 25). Just as in the Big Bang gravity evolves the 
initial density fluctuations imprinted on the cosmic microwave background (|26p 
into today's star and galaxy distributions fl27j), viscous hydrodynamics converts 
the initial energy density fluctuations of the "Little Bangs" created in heavy-ion 
collisions into a distribution of anisotropic harmonic flow components (I28ll241 l29 | 
30,31]). The low shear viscosity of the QGP allows traces of the initial fluctuations 
to survive to the final freeze-out stage, i.e. the flow anisotropics generated by 
them are not completely erased by dissipation and can be used to determine the 
QGP viscosity. While we have observational access only to a single instance of 
the Big Bang (the one that created the universe in which we live), RHIC and the 
Large Hadron Collider (LHC) have created billions of Little Bangs, each with a 
different distribution of initial fluctuations and its own hydrodynamic response 
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to it. From this extraordinarily rich data set both the QGP transport coefficients 
and the initial fluctuation spectrum can be reconstructed with precision. While 
this task is not yet complete, we review the dramatic progress and significant 
success achieved over the last few years. 

2 RELATIVISTIC HEAVY-ION COLLISION DYNAMICS 

2.1 Second order viscous relativistic fluid dynamics 

Hydrodynamics is an effective macroscopic description of a system that is in 
approximate local thermal equilibrium. It can be derived from the underlying 
microscopic (kinetic) description through an expansion in gradients of the local 
thermodynamic variables. In zeroth order of gradients one obtains ideal fluid 
dynamics. Navier-Stokes theory accounts for terms that are linear in gradients. 
These describe momentum and baryon number transport across fluid cells due to 
bulk and shear viscosity as well as heat conduction. Relativistic Navier-Stokes 
theory is, unfortunately, acausal (32): superluminal signal propagation causes 
short wavelength perturbations (which are outside the validity range of the gra- 
dient expansion) to grow exponentially, leading to physical and numerical insta- 
bilities. Any numerical implementation of relativistic dissipative fluid dynamics 
thus requires the inclusion of terms that are second order in gradients. Some 
of the associated transport coefficients have the physical meaning of microscopic 
relaxation times for the dissipative flows. Like the viscosities and the heat con- 
ductivity, they depend on the microscopic scattering dynamics; they encode the 
time delay between the appearance of thermodynamic gradients that drive the 
system out of local equilibrium and the associated build-up of dissipative flows 
in response to these gradients, thereby restoring causality. 

The most widely used theoretical frameworks for second-order dissipative rela- 
tivistic fluid dynamics are the 35-year-old Israel-Stewart theory (33) and its recent 
generalizations (j3^ l3l ^l3ini3Tl l3li | l3l ^l4T)j) . In its full glory, this theory contains a 
large number of second-order terms (|4ip , but none of the present numerical imple- 
mentations includes all of them. [Some exploratory studies including third-order 
gradient terms have also been published (jUJ, and they appear to improve the 
convergence of the gradient expansion.] While a full investigation is still outstand- 
ing, spot checks indicate that in practice not all second-order terms are equally 
important. Accounting for non-zero relaxation times at all stages of the evolution 
is, however, crucial: they limit excursions from local equilibrium, thereby both 
stabilizing the theory and improving its quantitative precision (|37y38p . 

For a system with small or zero net baryon number, heat conduction effects 
can be neglected. This simplification is generally being employed at LHC and 
top RHIC energies but will have to be abandoned at lower energies. This leaves 
the bulk and shear viscosities, £ and 77, and the relaxation times and t v for the 
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bulk and shear viscous pressures as the main transport coefficients controlling 
the collective dynamics. They can be usefully characterized by their unitless 
combinations with the entropy density s and the temperature T, resulting in 
the specific viscosities n/s and £/s and the scaled relaxation times Tt^^ v . For 
consistency of the theory the microscopic relaxation rates must be much larger 
than the scalar expansion rate 6 = d^u^ of the fluid, t„ c,0 <S 1. {vf'ix) is the flow 
4- velocity profile of the fluid.) 

2.2 Equation of state (EOS) 

Following recent progress in lattice gauge theory, the model equations of state of 
the past have been increasingly replaced by parametrizations of the lattice QCD 
EOS for which now nicely converged results are available (at least for baryon-free 
systems) (|4"3]) . Below T c ~ 155— 160 Me V (j4~3j) (the pseudo-critical temperature at 
which quarks and gluons combine into hadrons in a smooth but rapid crossover 
transition (|44p ) the lattice QCD EOS is matched to that of a hadron resonance 
gas (|45p . What matters for the development of collective flow is the relationship 
between the pressure p (whose gradients provide the accelerating forces) and 
the enthalpy e+p which embodies the inertia of the fluid (e is the comoving 
energy density). For ideal fluids (ignoring the viscous terms) the relativistic 
Euler equation can be written in the form ii u = , where we expressed the 

EOS p{e) approximately in terms of the squared speed of sound c 2 s = dp/de as 
p = c 2 s e. This shows that the evolution of flow in response to pressure gradients is 
controlled by the stiffness dp/de of the EOS, i.e. the key hydrodynamic ingredient 
is the speed of sound as a function of local energy density, c^(e). The connection 
between energy density e and temperature T (which depends on the number of 
active degrees of freedom in the medium and thus on its microscopic composition) 
is not needed for the hydrodynamic evolution; it is required, however, for the 
calculation of thermal electromagnetic radiation from the expanding fireball and 
final particle emission whose emission rates and energy distributions depends on 
the local temperature T. 

2.3 Quark-gluon plasma transport coefficients 

The calculation of QGP transport coefficients from QCD is difficult. Very little is 
known from first principles about heat conductivity in QCD matter. Leading or- 
der results in a weak-coupling expansion (<?<d) for the shear and bulk viscosities 
of a quark-gluon plasma were obtained in (|4"6ll4"T|) . For the coupling strengths 
relevant in a QGP at temperatures of a few T c , g ~ 1 — 2, these calculations 
provide not much more than an order of magnitude estimate. Taken literally, 
the leading-order perturbative value of the QGP shear viscosity (|46p is too large 
by a factor 3-10 to explain the experimentally measured large anisotropic flow 
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coefficients and the correspondingly required short QGP thermalization times 
(|48U49p . The extraction of r\ and £ directly from lattice QCD requires an analytic 
continuation from imaginary to real times. This is numerically costly and so far 
has only yielded results with large error bands (|50p . 

An alternate approach exploits the AdS/CFT correspondence relating strong- 
ly coupled conformal field theories (CFT) to classical gravity in weakly curved 
Anti-de-Sitter (AdS) space-time geometries (f5T]) . This approach has established 
approximate strong-coupling limits for the specific shear and bulk viscosities 
r?/s~l/(47r) d2H) andCA?^2(i-c;?) (J52J. For rj/s this is significantly below the 
perturbative QCD estimate while Q/s (which is proportional to the interaction- 
induced deviation from the conformal limit, \—c 2 s ) is much larger than in pertur- 



Below T c , a number of calculations of rj/s and £/s based on effective hadronic 
interaction models (see Kapusta (|49p for a review) and on transport models (|53p 
have been performed. Combining all calculations one finds that, generically, rj/s 
has a minimum near T c , rising steeply below and slowly above T c . The specific 
bulk viscosity £/s, on the other hand, peaks near T c (fMp . due to the breaking of 
conformal symmetry by long-range critical correlations; above and below T c , the 
ratio (,/rj generically decreases, allowing shear viscous effects to dominate over 
bulk viscous ones. The strongest effects from bulk viscosity (if any) are expected 
near T c . 

In view of limited theoretical knowledge of the QGP transport coefficients, 
recent attention has focussed on extracting them phenomenologically from ex- 
perimental data, by comparison with viscous fluid dynamical simulations. 

2.4 Hadron cascade stage, and chemical and kinetic freeze-out 

Below T c chemical reactions between different hadronic species become too slow 
to maintain chemical equilibrium in the exploding fireball (55,56). This leads to 
"chemical freeze-out" : While maintaining some degree of local kinetic equilibrium 
through quasi-elastic resonance scattering, the final stable hadron yields hardly 
change any more. An important exception are the baryon and antibaryon abun- 
dances which are somewhat depleted by annihilation in the hadronic rescattering 
stage (|57y58p . After correcting for baryon-antbaryon annihilation, the finally ob- 
served stable hadron yields reflect approximate chemical equilibrium abundances 
with a temperature T c h em close to T c (|58j) . 

If one continues with viscous hydrodynamics to describe the hadronic phase 
below T c , down to a final kinetic decoupling temperature T^ in < T c h em > one must 
account for chemical freeze-out at T c ^ em by assigning temperature-dependent non- 
equilibrium chemical potentials to each hadron species (113P59P - It was found that 
an incorrect chemical composition at Tkin can seriously distort the distribution of 
the hydrodynamically generated momentum anisotropy over the different hadron 
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species and over their transverse momenta f)13[) . Without accounting for the 
non-equilibrium hadronic chemical composition below T c hemi phenomenologically 
extracted values for the QGP shear viscosity can be wrong by 100% (|60y6ip . 

In spite of large resonant hadronic scattering cross sections, however, the col- 
lective expansion is so fast that even kinetic equilibrium is difficult to maintain 
below T c . Once in the hadronic phase, the fluid becomes so dissipative (|16p that 
the macroscopic hydrodynamic description becomes unreliable (|17p and should 
be replaced by a microscopic solution of the coupled Boltzmann equations for the 
various hadronic phase-space distributions ("hadron cascade", see e.g. (|62y63|) ). 
Still, many studies today continue to use the fluid dynamical short-cut all the way 
down to Tkj n . In such an approach, final hadron distributions are computed by 
converting the hydrodynamic output on a hypersurface of constant decoupling 
temperature Tk; n into hadron momentum distributions using the Cooper-Frye 
prescription, with non-equilibrium chemical potentials for the various hadron 
species and non-equilibrium corrections 5f(x,p) to the local phase-space distri- 
bution (|33|l64p that account for the viscous corrections to the energy momentum 
tensor on this decoupling surface. 

The more reliable, but much more expensive hybrid approach couples a viscous 
fluid dynamic description of the QGP phase with a microscopic Boltzmann simu- 
lation of the hadronic phase (|17U65p . VISHNU (|17p matches the (2+l)-dimensional, 
longitudinally boost-invariant viscous hydrodynamic algorithm VISH2+1 (|37U38p 
to the well-known UrQMD cascade ([62]), in MUSIC+UrQMD {(65|) the evolution is 
fully (3+l)-dimensional. In these hybrid approaches the Cooper-Frye algorithm 
is used to convert hydrodynamic output into particle phase-space distributions 
on a switching surface of constant temperature T c h e m) assuming chemical equi- 
librium yields (|17p . UrQMD then propagates these particles until all interactions 
cease and all unstable resonances have decayed. For good statistics, the final 
hadronic cascade stage is simulated many times for each hydrodynamic event. 

The influence of the dissipative hadronic stage on final observables increases 
relative to that of the early QGP stage as the collision energy decreases and 
the fireball spends less time as a QGP and a larger fraction of its history in 
the hadronic stage. Hybrid approaches such as VISHNU thus become more and 
more important as one moves down in energy. Some features measured at lower 
energies during the recent RHIC beam energy scan, such as different differential 
elliptic flow coefficients V2(pt) for baryons and antibaryons (|66p . are difficult to 
understand in pure hydrodynamics, and a hybrid approach may prove essential 
for their interpretation. 

At the highest available collision energies at the LHC, most of the anisotropic 
flow is created before hadronization, and the relative impact of the late hadronic 
stage is weaker (|67p . Good descriptions of the charged hadron p^-distributions 
and anisotropic flows v n and v n (pT) have been obtained both with pure viscous 



Flow and Viscosity in Relativistic Heavy Ion Collisions 



8 



hydrodynamics ([681 [691 [701 EI]) and with the hybrid code VISHNU (|7J) (see also 
recent work with MUSIC+UrQMD (65) and with hybrid codes coupling ideal fluid 
dynamics with UrQMD f)73|, IT4 [) ) . The hybrid approach yields a better overall 
description of the centrality-dependent balance between the evolution of radial 
and anisotropic flow, which particularly affects the £>T-spectra and differential 
elliptic flow of heavy particles such as protons (172p . 

2.5 Pre-equilibrium dynamics 

The hydrodynamic stage of ultra-relativistic heavy-ion collisions is preceded by a 
short (0.2— 1.5 fm/c) but very dense pre-equilibrium stage. During this stage the 
bulk of the energy density is contributed by gluon fields from the low-momentum 
components of the wave functions of the colliding nuclei. These fields are so 
intense that, in spite of weak QCD coupling a s = g 2 /{Ait) ~ 0.3 < 1, interaction 
mean free paths are too short for a quasiparticle-based kinetic approach to make 
much sense. A pre-equilibrium dynamical approach based on the non-Abelian 
evolution of classical gluon fields (the "Glasma" ([751 176 1 [77 1 178j) ) is probably a 
better starting point. 

Due to the finite size and anisotropic shape of the initial energy density distri- 
bution in the nuclear reaction zone, which fluctuates from collision to collision, 
any sort of interaction among the primordial QCD degrees of freedom in the 
fireball will cause non-vanishing radial and anisotropic flows even before the sys- 
tem has thermalized and viscous hydrodynamics becomes applicable. In fact, 
even if the system free-streams, correlations are generated between the average 
momenta and positions of the constituents that, when matched by the Landau 
prescription to a hydrodynamic form of the energy-momentum tensor, translate 
into non-vanishing anisotropic flow velocities on the matching surface. The hy- 
drodynamic stage thus starts with non-zero radial flow and non-vanishing flow 
anisotropics, but at the same time the pre-equilibrium evolution tends to dampen 
the initial spatial anisotropics and hence the anisotropic pressure gradients that 
drive anisotropic flow during the subsequent hydrodynamic stage. One therefore 
expects a certain degree of complementarity between pre- and post-equilibrium 
flow, resulting in reduced sensitivity of the final spectra to the starting time of 
the hydrodynamic stage (|711l79p . 

Model studies based on ideal fluid dynamics suggest that non-zero initial ra- 
dial flow at the beginning of the hydrodynamic evolution is phenomenologically 
preferred by some final state observables (|80p . However, faced with the con- 
ceptual difficulties of dealing with the early pre-equilibrium dynamics properly, 
most studies simply assumed an early hydrodynamic starting time with zero ini- 
tial transverse flow velocities. With the recent discovery of the importance of 
initial-state density fluctuations ("hot" and "cold" spots) in the early fireball 
this head-in-the-sand attitude became increasingly untenable. During the last 
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year, the dynamical IP-Glasma model (|81|) was developed which builds on the 
IP-Sat (Impact Parameter dependent Saturation) model (|82|) to generate finite 
deformed fluctuating initial gluon field configurations in the transverse plane, and 
then evolves them with classical Yang-Mills dynamics ([75l 176 1 177 1 178]) . While the 
lack of thermalization and of longitudinal fluctuations are still weaknesses of this 
model, it is the first semi-realistic approach to describing the pre-equilibrium 
stage dynamically, matching it consistently to the hydrodynamic staged Fig- 
ure [TJ shows three snapshots of the transverse energy density profile from this 
model. 



3 INITIAL-STATE DENSITY AND SHAPE FLUCTUATIONS 

3.1 Harmonic eccentricity and flow coefficients 

The development of anisotropic flow is controlled by the anisotropies in the pres- 
sure gradients which in turn depend on the shape and structure of the initial 
density profile. The latter can be characterized by a set of harmonic eccentricity 
coefficients e n and associated angles 3> n : 

j$ x _ f r dr dcj)r 3 e i 't' e(r,<fi) m$„ _ f r dr dip^e 1 ™^ e(r,4>) ^ 
j r dr d<pr 3 e(r,(j)) ' n j r dr d(j)r n e(r,(j)) 

(1) 

where e(r, <j)) is the initial energy density distribution in the plane transverse to 
the beam direction. When, for collisions between nuclei of the same species, e is 
averaged over many events and the angle eft is measured relative to the impact 
parameter vector, there is a symmetry between (j) and —0 as well as between <j> 
and 4> + tt, and all odd e n coefficients vanish. 

An important insight (123 1 124 11251) has been that, due to event-by-event fluctu- 
ations of the transverse positions of the nucleons inside the colliding nuclei ()22j) . 
and of the gluon density profiles inside those nucleons ([75|I76| [771 ITHjIHTj IM1 I85 ]) 
(see Figure [TJ), these symmetries do not hold in an individual collision event. 
Therefore, in every collision all eccentricity coefficients are usually non-zero, driv- 
ing anisotropic flow components of any harmonic order whose magnitudes and 
directions fluctuate from event to event. The statistical distributions of e n and 
<3? n which, in a hydrodynamic picture, control the statistical distributions of the 
final anisotropic flows v n and their directions Vl/ n , are of quantum mechanical 
origin and depend on the internal structure of the colliding nuclei (see Sec. I3.3|) . 

The anisotropic flow coefficients v n and their associated flow angles ^f n are 

1 It has been suggested that, since classical Yang-Mills dynamics does not lead to local ther- 
malization, a different matching scheme (83) should be used that, unlike Landau matching, does 
not rely on small deviations from local equilibrium. This interesting suggestion still needs to be 
fully worked out for fluctuating initial conditions. 
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defined in analogy to Eq. ([1]) as 

v n e A = (e m ^), v n ( PT )e in *^ = (e ini ?) PT , (2) 

where ip is the azimuthal angle of the transverse momentum vector pt of the 
emitted particle in the laboratory frame. Note that the v n characterize a sin- 
gle collision. In the left equation, which defines the "pr-integrated" or "total" 
anisotropic flow, the average is over all particles of a given kind (identified by 
species or not) in the event. The right equation defines the "pr-differential" flow 
and averages only over particles with a given transverse momentum. v\, V2, and 
^3 are known as directed, elliptic, and triangular flow, respectively. 

Due to the limited number of particles emitted from a single event, statis- 
tically precise measurements of v n and ty n can only be obtained from particle 
distributions that have been averaged over many events. One defines the mean 
flow coefficient v n and the mean flow angle ^ n through a Fourier decomposition 
of the experimentally determined, event averaged particle distribution: 

= ^ f 1 + 2 Xj *n cos(n(,-*„)) , (3) 

where N = (N) is the mean number of particles of interest per event (charged 
hadrons or identified particles of a specific species). The generalization to pr- 
differential v n (pT) and ^f n {pr) in terms of dN / {pxdpTdip) is obvious. The Fourier 
coefficients are given by: 

v n = ((cos[n(^-*„)]», or equivalents v n = ((e^))^*™, (4) 

where ((. . .)) denotes the double average over all particles in an event and over 
all events. Note that, due to event-by-event fluctuations of the flow angle *$> n , 
the right expression is not identical with the average of Eq. ([2]) over collision 
events as is usually assumed. Observable consequences of these event-by-event 
fluctuations of the flow planes for generally employed experimental v n measures 
(to be discussed next) are only now being investigated (|86 H87j) . 

Since the flow planes are not experimentally known, the anisotropic flow coef- 
ficients are calculated using azimuthal angular correlations between the observed 
particles. In the case of two particle correlations l^e in ^ tpi ~ <( "^J) the measurement 
is proportional to (v%)' under the assumption that the only azimuthal correlation 
between particles is due to the common correlation with the flow plane this corre- 
lator can be factorized into (( e in ( < Pi-*»))( e - in ^ 2 -*™))) = (v%). From this method, 
the experimentally reported anisotropic flow coefficients are therefore obtained 
as the root mean square value yj (v'fy. Due to event-by-event fluctuations in the 
anisotropic flow, the event averaged / {v n ) k for k > 2. We denote by 
v n {2} = \J (v'n) the anisotropic flow extracted from two-particle correlations. 

In practice, not all azimuthal correlations in the data are of collective origin. 
Additional "non-flow" correlations arise from resonance decays, jet fragmenta- 
tion, and Bose-Einstein correlations ()88|) . They can be suppressed by appropriate 
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kinematic cuts (|88p or by using multi-particle correlations known as higher-order 
cumulants ([89]) . v n {4} and v n {6} denote the anisotropic flow coefficients obtained 
from the fourth and sixth order cumulants, respectively. 

In addition to being less sensitive to non-flow contributions, the higher order 
cumulants, which involve higher moments of the event-by-event v n distribution, 
also depend differently on the variance cr Vn of that distribution. If a Vn <C (v n ), 
one finds up to order (a Vn / {v n }) 2 




This illustrates that the difference between v n {2} and v n {4} is sensitive to the 
width of the v n distribution, in addition to non-flow effects. After correcting 
v n {2} for non-flow (|90p . this difference can be used to estimate the variance of 
the event-by-event flow fluctuations (see Sec. 16. 4Q . 

3.2 Centrality classes 

Heavy ions are extended objects, and the system created in a head-on collision 
is different from that in a peripheral collision. We therefore categorize nuclear 
collisions by their centrality, parametrized by the impact parameter b which is, 
however, not a direct observable. Experimentally, the collision centrality can be 
inferred from the number of produced hadrons, if one assumes that this multi- 
plicity is a monotonic function of b. Knowing which fraction of the total hadronic 
cross section is observed in the experiment one can divide the measured distri- 
bution of produced particles in centrally classes corresponding to the percentile 
of total hadronic cross section. 

In addition to impact parameter or fraction of total hadronic cross section, one 
can also use the so-called number of wounded nucleons or the number of binary 
nucleon-nucleon collisions to characterize the collision centrality. These quantities 
are defined within the Glauber model (|9ip and are stochastic functions of the 
impact parameter. Phenomenologically it is found that soft particle production 
scales roughly with the number of participating nucleons whereas hard processes 
scale with the number of binary collisions. 

Anisotropic flow is not measured in a single event but in a centrality class. 
Therefore, event-by-event fluctuations due to impact parameter fluctuations within 
a centrality class will add to the initial-state fluctuations mentioned in Sec. 13.11 
(and further discussed in the following subsection) to determine the spectrum of 
final-state flow (v n and fluctuations. 

3.3 Models for the primordial fluctuation power spectrum 

Various theoretical approaches have been used to model the initial energy and 
entropy density profiles. The most common approach is the Monte Carlo (MC) 





and 
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Glauber model (|9ip . In this model the positions of the nucleons inside the two 
colliding nuclei are sampled according to the measured nuclear density distribu- 
tion, accounting for their finite size (92). The sampling procedure introduces 
event-by-event fluctuations of the nucleon positions, representing the quantum 
mechanical fluctuations of the outcome of a position measurement on the nu- 
cleons whose probability distribution is given by the (smooth) nuclear ground 
state wave function. The nucleons travel on straight-line trajectories and collide 
if their distance in the transverse plane is smaller than the radius corresponding 
to the total inelastic nucleon-nucleon cross-section. It is assumed that this radius 
is independent of the number of interactions the nucleons already had. Nucle- 
ons in the target and projectile that have had at least one interaction are called 
participants or wounded nucleons. Soft particle production is assumed to be 
proportional to the number density of wounded nucleons whereas hard (high-p^) 
processes are taken to scale with the number density of binary nucleon-nucleon 
collisions. The initial entropy or energy density profile is typically taken pro- 
portional to a linear combination of the wounded nucleon and binary collision 
densities which are strongly fluctuating in the transverse plane and from event 
to event. 

In the Monte Carlo Kharzeev-Levin-Nardi (MC-KLN) model (|93|) the entropy 
production is determined by the initial gluon production, calculated by pertur- 
bative merging of two gluons from the projectile and target nuclei where (in the 
spirit of Color Glass Condensate ideas (i94|) ) the (pr-unintegrated) gluon struc- 
ture functions are parametrized by a position dependent gluon saturation mo- 
mentum Q s ([95]) . Q s (r±) is computed from the longitudinally projected density 
of wounded nucleons whose positions are sampled as in the MC-Glauber model. 

The MC-Glauber and MC-KLN models do not account for fluctuations of the 
gluon fields inside the colliding nucleons which are characterized by a correlation 
length ~ 1/Qs i n the transverse plane (|84p . For this reason, they cannot re- 
produce the measured multiplicity fluctuations in pp collisions (I85p . Gluon field 
fluctuations can be imprinted on the MC-Glauber or MC-KLN profiles a poste- 
riori using an algorithm developed in (|96p . In the IP-Glasma model, described 
at the end of Sec. 12.51 gluon field fluctuations are imprinted on the IP-Sat model 
(82) using ideas from the Color Glass Condensate/Glasma theory fj75 |, l76 ^l94[) and 
then evolved using classical Yang-Mills dynamics (177}l78p to a matching surface 
after which hydrodynamics takes over. 

DIPSY is a Monte Carlo event generator based on gluon radiation from colored 
dipoles (via dipole splitting) using BFKL evolution ([97]) . The nucleons in the 
colliding nuclei are described by a triangle of color dipoles whose position is 
sampled from a Woods-Saxon distribution just as in the MC-Glauber and MC- 
KLN models. Gluon density fluctuations arise from the Monte Carlo sampling of 
the radiative gluon shower. 
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The Monte Carlo approach used in these models allows for an event-by-event 
calculation of the eccentricities e n and thus for the determination of the higher 
order moments of the e n distributions. Because of the non-zero widths of these 
distributions, the measured final flow power spectrum obtained from a given 
moment of the v n distribution, say v n {k}, should be compared with the initial 
fluctuation power spectrum e n {k} obtained from the analogous moment of the 
e n distribution (|88|) . 

Figure [2] shows the (e n ) power spectrum of the initial energy density for the 
MC-Glauber, MC-KLN and IP-Glasma models, for four centrality classes. In 
the most central collisions, — 0.2%, the e n are entirely due to fluctuations, 
and therefore all (e n ) have roughly equal magnitudes. Proceeding to less central 
collisions, the nuclear overlap region develops a pronounced elliptic geometric 
deformation which increases (£2) (and, to some extent, also (£4) and (£6)) much 
more strongly than the odd eccentricity coefficients which remain fluctuation 
dominated. Importantly, Figure [2] reveals a significant model dependence of 
the {e n ) spectrum. This implies on the one hand that the extraction of QGP 
transport coefficients (see Sec. 14. 2|) is complicated by uncertainties related to 
initial model ambiguities, but on the other hand that, by constraining the theory 
with a sufficiently large set of anisotropic flow and correlation measurements (see 
Chap. [6]), one can discriminate phenomenologically between different initial state 
models and identify the one that produces the correct initial fluctuation power 
spectrum. 

The different initial state models also produce different correlations between 
the "participant plane angles" $ n associated with e n . For the MC-Glauber and 
MC-KLN models these were studied in (|98l [99l HOOD . The evolution of these 
initial-state participant plane correlations through non-linear viscous hydrody- 
namics into final-state flow-plane correlations (jlOOtllOi]) and the experimental 
measurement of the latter (|102yi03|l will be discussed in Sec. 16.31 

4 HYDRODYNAMIC RESPONSE TO INITIAL-STATE FLUC- 
TUATIONS 

4.1 Viscous effects on radial and anisotropic flow 

The viscosity of a fluid is related to its ability to return to local thermal equi- 
librium after being driven away from equilibrium by gradients in its macroscopic 
flow pattern. When the corresponding relaxation times approach zero, viscous 
effects disappear. Small viscosities are therefore related to short relaxation times, 
i.e. strong interactions among the microscopic constituents. 

For non-zero relaxation times, viscous effects cause the microscopic momen- 
tum distribution in the local fluid rest frame to deviate from its exponential local 
equilibrium form. Bulk viscosity causes locally isotropic deviations from equi- 
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librium, adding a diagonal contribution, the bulk viscous pressure lit?- 7 , to the 
stress tensor T %3 in the local rest frame. In Navier-Stokes approximation (which 
ignores retardation) it is proportional to the scalar fluid expansion rate 6 at the 
location of the fluid cell, H(x) = —Q9(x). For an expanding fluid this is negative, 
i.e. the bulk viscous pressure counteracts the expansion. In an isotropically ex- 
panding fireball, bulk viscosity reduces the radial acceleration and thus inhibits 
the buildup of radial flow. 

Shear viscosity causes locally anisotropic deviations from equilibrium, resulting 
in an anisotropic contribution to the local rest frame stress tensor. This shear 
viscous pressure is driven by shear flow and acts against flow anisotropics. In 
relativistic heavy-ion collisions, due to approximate boost-invariance along the 
beam direction of the physical processes that generate the quark-gluon plasma, 
the initial expansion rate is highly anisotropic and much larger along the beam 
direction than transverse to it. Transverse flow builds only later, in response to 
transverse pressure gradients in the initial state. The main effect of shear vis- 
cosity is that it tries to equalize the expansion rates along different directions, 
by building up a shear viscous pressure tensor ir 13 in the local rest frame that 
reduces the longitudinal and increases the transverse pressure. Less work is done 
by longitudinal pressure, causing the QGP to cool less rapidly (at least initially 
when cooling is dominated by longitudinal expansion) but simultaneously in- 
creasing the build-up of transverse flow. Transverse anisotropics of the initial 
fireball geometry, reflected in anisotropic transverse pressure gradients, generate 
anisotropics in the developing transverse flow; shear viscosity reduces these flow 
anisotropics, i.e. it degrades the medium's ability to convert initial transverse 
pressure anisotropics into final transverse flow anisotropics. 

Anisotropic transverse flow influences the shape of the transverse momentum 
spectra of the finally emitted particles through a direction-dependent blue-shift 
factor. By reducing the azimuthally symmetric component of transverse flow, 
a.k.a. radial flow, bulk viscosity leads to steeper p^-spectra; by boosting radial 
flow, shear viscosity renders them flatter. Numerically, it was found (|104p that 
in heavy-ion collisions shear viscous effects dominate over bulk viscous ones by 
about a factor 5-10. Experimentally, they are not easy to separate since both 
types of effects modify the slopes of the p^-spectra as well as the anisotropic flow 
coefficients discussed below. A clear strategy for systematically isolating shear 
from bulk viscous effects still needs to be formulated. Lacking such a strategy, 
most researchers presently allow themselves to be guided by available theoretical 
studies and interpret anisotropic flow measurements entirely in terms of shear 
viscosity, i.e. they ignore bulk viscosity in their theoretical models. 
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4.2 Extracting the QGP shear viscosity from experimental data 



Shear viscosity affects the p-p-spectra of the finally emitted particles in two dis- 
tinct ways: (i) it increases the magnitude and decreases the anisotropics of the 
hydrodynamically generated transverse flow, and (ii) it causes a deviation Sf of 
the final phase-space distribution f(x,p) from its isotropic local equilibrium form 
fo(x,p): f(x,p) = fo{x,p) + 5f(x,p). The first of these two effects reflects the 
amount of shear viscosity over the entire expansion history of the fireball; the 
second effect, 5f, is only sensitive to the shear viscosity at the final conversion to 
hadrons (assuming, as required for the applicability of viscous hydrodynamics, 
that microscopic relaxation rates are much larger than the macroscopic expan- 
sion rate) . 5f increases with transverse momentum roughly as p T where a ranges 
between 1 and 2 (JM])- 5f is thus small at low p T <5T (where T ~ 100-160MeV 
is the decoupling temperature), but at larger p^>2 — 3GeV it becomes so big 
(|37|) that the near-equilibrium expansion (and thus the hydrodynamic prediction 
for the p^-distribution) can no longer be trusted. Being predominantly interested 
in the shear viscosity (ji/s)qqp of the early QGP phase, one would like to mini- 
mize in the analysis the effects from Sf which reflect only the late hadronic shear 
viscosity. This can be achieved by studying pr-integrated observables which are 
dominated by transverse momenta py < 5 T where Sf is small. 

A second consideration is that the hydrodynamically generated transverse mo- 
mentum anisotropy is distributed in the measured final state over a large number 
of different hadronic species. Since heavier particles are affected by radial flow 
more strongly than lighter particles (1105[) . resulting in flatter py-spectra at low- 
px where most final hadrons find themselves, the chemical composition of the 
final state, in particular the light-to-heavy particle ratios, control where in px 
the hydrodynamic momentum anisotropy ends up: For heavy hadrons, radial 
flow pushes the flow anisotropies to larger px; at low px the anisotropic flow is 
dominated by light hadron species. Extracting the shear viscosity rj/s from the 
pr-dependent anisotropic flow coefficients thus introduces an undesirable fragility 
related to the precision with which the theoretical model describes the chemical 
composition of the fireball at freeze-out and the shapes of the different hadronic 
pr-spectra that control the distribution in px of the hydrodynamically gener- 
ated momentum anisotropy. A more robust extraction of r/ / s uses the azimuthal 
anisotropy coefficients of the pr-integrated angular distributions, summed over 
all hadronic species (|lU6yiU7lllU8yiU9|) (the so-called "charged hadron i> n "), and 
later checks the p^-differential flow coefficients v n (px), for all charged hadrons 
together and for individual hadron species separately, as additional tests whether 
the model also correctly describes the relative hadron abundances and their px- 
distributions in the final freeze-out stage. 

Figure [3] shows an attempt of such an extraction of the specific QGP shear 
viscosity from charged hadron elliptic flow (t^) data collected at RHIC. The plot 
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shows v 2 scaled by the initial fireball ellipticity e = Ei (e characterizes the elliptical 
spatial deformation of the initial pressure gradients that drive the hydrodynamic 
expansion of the fireball) as a function of collision centrality, represented by the 
charged hadron multiplicity density per unit overlap area (l/S)(dN c \ l /dy) (the 
most central collisions have the largest multiplicity density). The theoretical 
curves were obtained with the VISHNU hybrid model (117h with a temperature- 
independent specific QGP shear viscosity (j?/s)qgp- The model results show a 
monotonic dependence on (t]/s)qqp and nicely reproduce the shape of the mea- 
sured centrality dependence. By normalizing V2 by £2, the theoretical curves 
become insensitive to the experimental method used for measuring V2 as long as 
the same procedure used in determining vi is also used to compute the initial 
ellipticity £2 (|108U109]) . However, the magnitude of £2 differs by about 20% be- 
tween the two initial-state models studied in the left and right panels of Figure [31 
shifting the normalized experimental data up or down relative to the theoretical 
curves by similar amounts. As a result, the (??/s)qgp value extracted from the 
comparison is uncertain by a factor 2—2.5 and depends on the specific choice of 
the initialization model. In Sec. we review how this model ambiguity can be 
resolved by using experimental information from all flow harmonics v n and their 
systematic dependences on collision energy and centrality and on system size. 

5 TRANSVERSE MOMENTUM SPECTRA AND RADIAL FLOW 

5.1 Radial flow systematics at RHIC and LHC energies 

Radial flow causes a blueshift of the transverse momenta of the finally emitted 
hadrons, leading to flatter px and mx distributions, especially at low px where 
non-relativistic kinematics allows to express px ~ Px~ + fn(vx) as the sum of 
a thermal contribution p^ (which depends only on the decoupling temperature 
7km and is independent of the hadron mass m) and a collective flow component 
with average flow velocity {vt) which is proportional to the hadron mass (|105p . 
This is shown in Figure [4j Due to collective flow, the heavy-ion collision spectra 
in the right panel are flatter than the pp spectra in the left panel, and due to 
~ 10% stronger flow at the higher collision energy, they are flatter at the LHC 
than at RHIC. For heavy hadrons such as protons, the radial flow generates a 
"shoulder" in the spectrum at low px which is more pronounced at LHC than at 
RHIC energies. In contrast, the soft parts (px < 2.5 GeV/c) of the pp spectra do 
not feature the flow-induced splitting between kaon and proton slope parameters 
that characterizes the heavy-ion collision spectra^ 

2 The different slopes of pion and kaon spectra in pp collisions have a different origin, unre- 
lated to flow: a large fraction of the pions in the final state arises from the decay of unstable 
heavier resonances, and these decay pions accumulate preferentially at low pr, due to kinematic 

constraints (1116). The larger inverse slope of the pp collision spectra at the higher LHC energy 
is an initial state effect, arising (among other factors) from a higher gluon density probed at 
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Blast-wave model fits (|105p of the hadron spectra in terms of two parameters, 
the decoupling temperature T^ in and the average transverse flow (vx) at kinetic 
decoupling, show that (vx) increases while T^m decreases with growing collision 
energy and as the collisions become more central ( f|118|lll9p and personal commu- 
nication by L. Kumar). This is consistent with expectations from a kinetic theory 
description of the decoupling process (j55j). Hydrodynamic simulations show that 
the stronger radial flow at higher energies and in more central collisions is mostly 
a consequence of higher initial fireball densities, leading to a longer fireball life- 
time. Stronger flow increases the mean px of the emitted hadrons while a lower 
decoupling temperature decreases it. Figure [5] shows that, at fixed collision 
energy, (px) grows as the collisions become more central; a similar analysis for 
central collisions as a function of collision energy shows that (px) also grows 
with collision energy. The positive effect on (px) from increasing radial flow thus 
dominates over the negative effect from the accompanying decrease of Tki n . 

The mean px of charged hadrons increases with increasing shear viscosity but 
decreases with increasing bulk viscosity. For Glauber model initial conditions it is 
also affected by the width w of the Gaussian smearing profile (|12U|ll2ip : increasing 
w produces smoother initial density profiles and reduces the radial flow and thus 
(pt)- Lowering the decoupling temperature Tkj n increases the fireball lifetime 
and thus the radial flow and (px) (|120p . 

Figures [3] and [B] show px spectra for pions, kaons and protons from Au+Au 
and Pb+Pb collisions at RHIC and LHC. As seen in the right panel in Fig- 
ure |H hydrodynamic models with (HKM (JT^D and MUSIC+UrQMD ([65])) and with- 
out a hadronic cascade afterburner (VISH2+1 (j68[) and Krakow (|70p ) describe the 
data quite well; hydrodynamic models without hadronic afterburner that imple- 
ment chemical freeze-out at T c h em ~ T c overpredict the proton yields, by ignoring 
baryon-antibaryon annihilation in the hadronic phase, but still reproduce the 
shape of the spectra. As seen in the left panels of Figure [6j the azimuthally 
averaged px spectra cannot distinguish by themselves between different initial 
conditions and different values of the QGP shear viscosity: While larger 77/ 's val- 
ues cause stronger radial flow, this can be compensated for by assuming longer 
thermalization times and a corresponding later start of the hydrodynamic evo- 
lution stage (108) which shortens the time available for building radial flow. 
A determination of the QGP transport coefficients thus requires the simulta- 
neous investigation of the azimuthally averaged p^-spectra together with their 
azimuthal anisotropics (see Sec.[6|). 

5.2 Radial flow fluctuations 

Initial-state fluctuations affect not only the shape but also the size of the ini- 
tial fireball (|120p . Figure \7\ shows initial wounded nucleon distributions for two 

LHC energies compared to RHIC (|117[> . 
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Au+Au collisions with equal numbers of wounded nucleons (N w = 100), but 
rms radii that differ by more than 40% (3.14 and 2.38 fm, respectively (|120p ). 
The larger pressure gradients associated with the more compact configuration 
lead, through hydrodynamic evolution, to a 10% larger mean pt of the finally 
observed charged hadrons (622 vs. 563MeV/c, respectively). The normalized 
variance c r /(r) of the fluctuations in the initial fireball radius increases from 
about 2—3% in central Au+Au collisions to more than 15% in peripheral colli- 
sions (|120p . Event-by-event viscous hydrodynamic evolution of these fluctuating 
initial states allows to compute the covariance of the resulting final state px fluc- 
tuations. For Glauber initial conditions evolved with rj/s = 0.08 and a hadronic 
bulk viscosity of Q/ s = 0.04, the theoretical predictions (|120p compare well with 
available experimental data from the PHENIX and STAR Collaborations (|122[) . 
While radial flow fluctuations are perhaps not the best observable for a precision 
measurement of the QGP viscosity, they provide a valuable consistency check for 
the hydrodynamic approach and can help to constrain the spectrum of initial- 
state fluctuations. 

6 ELLIPTIC AND OTHER ANISOTROPIC FLOW COEFFI- 
CIENTS 

6.1 Elliptic flow systematics 

Experimentally, the most direct evidence of hydrodynamic behavior comes from 
the observation of anisotropic flow ()123j) . The largest of the anisotropic flow 
coefficients is t>2, the elliptic flow. Like all v n , it depends on rj/s: Larger shear 
viscosity quickly reduces them. Therefore, the large elliptic flow observed at 
RHIC energies provides compelling evidence for strongly interacting matter that 
behaves like an almost perfect liquid (|118pi24|) . 

The viscous effects reducing the magnitude of the elliptic flow depend on the 
size of the system (38) which again depends on the collision centrality. The 
centrality dependence of the elliptic flow thus is an observable sensitive to the 
magnitude of rj/s. In Figure [2] (Sec. I4.2p we showed that the centrality depen- 
dence of ^2/^2 is nicely described by viscous hydrodynamic calculations. However 
the magnitude of rj/s used in these calculations should be considered as an av- 
erage over the temperature history of the expanding fireball since we know from 
other fluids that rj/s depends on temperature. In addition, we also know that 
part of the elliptic flow originates from the hadronic phase. Therefore, knowl- 
edge of both the temperature dependence and the relative contributions from the 
partonic and hadronic phases is required to quantify (t//s)qgp of the partonic 
fluid. 

At RHIC energies, the dissipative hadronic phase significantly affects the finally 
observed elliptic flow, complicating the accurate determination of (t]/s)qqp. In 
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Pb+Pb collisions at the LHC the higher collision energies produce a system that 
is hotter and has a longer-lived partonic phase. Consequently the hadronic con- 
tribution to the elliptic flow decreases, and this reduces the uncertainty on the 
determination of (t?/s)qgp- Measurements of the energy dependence of elliptic 
flow at both RHIC and at the LHC allow to vary in a systematic manner the 
contribution from both phases and probe the temperature dependence of rj/s 

Figure [8k shows the measured integrated elliptic flow at the LHC in one cen- 
trality bin, compared to results from lower energies. It shows that there is a 
continuous increase in the elliptic flow from RHIC to LHC energies. In compari- 
son to the elliptic flow measurements in Au+Au collisions at ^/s NN = 200 GeV, V2 
increases by about 30% at y/s^ = 2.76 TeV. Hydrodynamic models (1127012801291) 
and hybrid models (I16 pi30p that successfully describe flow at RHIC indeed pre- 
dicted an increase of ~10-30% in vi at the LHC. 

At low pr, not only the i>2 of charged particles but also that of identified 
particles at RHIC and LHC is well described by viscous hydrodynamics. The blue 
shift of the transverse momenta which depends on the particle mass also generates 
the characteristic mass splitting observed in a plot of V2 versus px for particles 
of different mass (|118tll31j) . Figure [8J3 shows the pion and proton elliptic flow 
measured by STAR compared to VISHNU model calculations (|109p . In the initial 
viscous hydrodynamic stage rj/s is taken to be temperature independent. The 
rj/s values required to describe the pion and proton p^-differential elliptic flow 
data are found to be one or two times the KSS bound (|2Tj) for Glauber or KLN 
eccentricities, respectively, in agreement with the values required to describe the 
PT-integrated total charged particle v-i- 

6.2 Triangular flow and higher flow harmonics 

The produced system responds as a fluid to the pressure gradients and converts 
the coordinate space distributions, characterized by e n , to long range momentum 
space correlations between the produced particles. In the last years it was realized 
that novel long range correlations first observed at RHIC, known as the "ridge" 
and "Mach cone", are, in fact, manifestations of anisotropic flow (|24 ^ll32tllQ2|) . 

Because the created matter distribution in a collision is inhomogeneous, n cov- 
ers a large range (Figure [2]). However, the shear viscosity reduces differences 
between the expansion velocities and therefore dampens the anisotropic flow co- 
efficients v n , and more strongly so for larger n. As a consequence, the magnitude 
and transverse momentum dependence of the v n coefficients provide a large set 
of observables (I102lll33pi34lll35lll36p to check the hydrodynamical paradigm and 
are, within this description, very sensitive to the magnitude of r\j s (|137p69p . 

In Figure [9k the measured V2, v$ and V4 are shown as functions of cen- 
trality (|7ip . These form the flow power spectrum created by hydrodynamics 
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in response to the initial fluctuation power spectrum shown in Figure [2j pj>- 
integrated charged hadron flow coefficients for n = l and n > 4 are not yet avail- 
able experimentally but have been predicted theoretically (|7ip and are expected 
to complement the Little Bang power spectrum in the near future. Clearly v$ 
and U4 have a rather weak centrality dependence while t>2 changes rapidly. This 
mirrors the spectrum of e n in Figure [2] which also shows the strongest central- 
ity dependence for n = 2. Within hydrodynamics, the initial e n power spectrum 
together with the shear viscosity (rj/s)(T) completely determine the centrality 
and pt dependences of the anisotropic flow coefficients v n , as well as their event- 
by-event fluctuations. The IP-Glasma initial conditions (|8ip . together with an 
average value for r]/s of 0.2 for Pb+Pb collisions at the LHC and a somewhat 
smaller value of 0.12 for Au+Au collisions at top RHIC energies, provide a good 
description of all presently available data for charged hadron v n and f n (pr) 
Examples are shown in Figure [9j As discussed in (j71[) , the measurements at both 
collision energies are also compatible with a temperature dependent specific shear 
viscosity (rj/s)(T) (|125p that has a minimum value of ^ = 0.08 at T c and rises 
moderately above and more steeply below T c (c./. the two curves shown in com- 
parison with V2{pt) data in Figure [HJa). On the other hand, initial fluctuation 
power spectra obtained from the MC-Glauber and MC-KLN models (also shown 
in Figure [2]) can not reproduce the measured flow power spectrum (C. Shen 
and Z. Qiu, private communication, and related discussion in (|138p ). Although a 
good simultaneous description of charged hadron V2 and v 3 can be obtained with 
rj/s = 0.08 and MC-Glauber initial conditions (fT39j) (but not with MC-KLN ini- 
tial conditions), both models fail when tested on higher flow harmonics and the 
widths of their event-by-event fluctuation distributions. 

Fluctuations also induce a non-zero directed flow v\ at midrapidity (I14ipi40"j) . 
It is strongly constrained by global momentum conservation (1140]) (which forces 
the pT-weighted directed flow to vanish at midrapidity) and has been measured 
both at RHIC (fT42j) and LHC (fT45j) . Its effect is clearly visible in Figure [Has 
a dipole shift of the matter density. Once still existing differences between the 
measurements have been clarified, this observable will form the bottom end of 
the flow power spectrum of the Little Bang. 

6.3 Flow angle correlations 

The anisotropic flow coefficients V2 and V3 are to a good approximation linearly 
proportional to the £2 and £3, respectively (|3ip. However, due to the nonlin- 
earity of the hydrodynamic evolution equations, the hydrodynamic response to 
the initial-state fluctuations spectrum is characterized by mode-mixing between 
different order flow harmonics which becomes large in mid-central and peripheral 
collisions where the large geometric ellipticity £2 drives a large elliptic flow V2 
which mixes with the other flow harmonics (|3iyi44pl0ip . Anisotropic flow coeffi- 
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cients such as and v§ thus depend on the magnitude of V2 and V3, respectively, 
and for mid-central collisions the flow angles ^4 and ^5 are uncorrelated with the 
participant plane angles $4 and $5 associated with £4 and £5 (f3"T|) . This charac- 
teristic change between the initial coordinate space e n and final momentum space 
v n power spectra and their associated angles is a strong test of the hydrodynamic 
paradigm and provides additional constraints on the specific shear viscosity and 
initial density fluctuation spectrum (jlOOp . 

Correlations between the flow planes ^ n can be measured either by using com- 
binations of two particle correlations to estimate each ^f n , or by using multi- 
particle cumulants. Figure 110a shows the ALICE measurements (jT02) of 5- 
particle cumulants that are sensitive to the correlation among different order 
flow planes. The correlator ((cos[3(^i+<P2) — 2(</?3+¥>4+¥>5)])) (where ((...)) de- 
notes the double average over the particles 1-5 within an event and over all 
events) is sensitive to (v^v^ cos[6( , 3 / 3— ^2)]) (where (. . .) denotes an average over 
events); the data show that ^3 and ^2 are uncorrelated for central to mid-central 
collisions. However, the three-plane correlation between the angles VPi, ^2 an d 
\p3, obtained by measuring the five-particle cumulant shown by the blue mark- 
ers, is already significant for mid-central collisions, in qualitative agreement with 
expectations from (linearized) hydrodynamic response (|14ip . 

In general, however, linearized hydrodynamic response is not sufficient. In Fig- 
ure 110b we show several correlation functions between flow angles corresponding 
to different harmonics as measured by ATLAS (|103p . plotted versus centrality, 
with central collisions (large -/Vp art ) on the the right and peripheral collisions 
(small Ap art ) on the left. The comparison to viscous hydrodynamic calculations 
(jlOOp with initial energy density profiles from the MC-Glauber and MC-KLN 
models (solid and dashed lines) shows good qualitative overall agreement; the 
corresponding correlations between the eccentricity planes <J? n in the initial state, 
on the other hand, behave quite differently - in magnitude, in their qualitative 
centrality dependence, and even in sign (jlUUp . The final-state flow angle correla- 
tions thus cannot be understood, even qualitatively, without taking into account 
the nonlinear hydrodynamic response of the QGP fluid to the fluctuating initial 
states (jlOOyiOlj) . The strength of these nonlinear effects and, in some cases, even 
the shape of the centrality dependence of these flow angle correlations were shown 
to be sensitive both to the details of the initial-state fluctuation spectrum and to 
the specific shear viscosity rj/s (jlOOp . 

6.4 Eccentricity and flow fluctuations 

As explained in Sec. 13. 1} experimentally the anisotropic flow coefficients are esti- 
mated from measured angular correlations between emitted particles. Different 
correlation functions probe different moments of the v n distributions, but none of 
them measures directly the mean v n . As seen in Eq. ([5]), the mean and variance 
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of v n can be approximately obtained from the two- and four-particle cumulants 
(after correcting the former for non-flow effects (f90|) ) as v n ph \J {v'f l {2}+v'^ l {A}) /2 
and o Vn ~ \J (v'n{2}— t>^{4}) /2jf| The blue markers in Figure II lb show the re- 
sulting experimental estimate for the normalized elliptic flow variance cr V2 /v2- It 
is found to be large (~40%), hence the shown V2 and only rough esti- 

mates. Since, for not too large impact parameters (|3ip . the magnitude of the 
elliptic flow is proportional to the eccentricity £2 of the initial nuclear overlap re- 
gion, we expect (and the top panel in Figur eTTb confirms) that the normalized 
variance of V2 should be close to that of £2, shown by the lines in Figure 1111 
for MC-Glauber (jHH) and MC-KLN (|95]) initial conditions. The difference be- 
tween the dash-dotted and solid lines, showing the exact ratio a e2 /£~2 from the 
respective model and its small-variance approximation, indicates the quality of 
the approximation used in the experimental extraction of a V2 / V2 ■ The rise of the 
experimental estimate for a V2 /v2 in central collisions is seen to be mostly due 
to this approximation. Although the centrality dependences of a £2 /£2 from the 
models and of cr V2 /v2 from experiment show clear similarities, neither model is 
able to match the data over the full centrality range. 

For a more detailed comparison the full event-by-event distribution of v n can 
be reconstructed experimentally, by removing as much as possible non-flow cor- 
relations and then unfolding the measured correlations with a response function 
that accounts for the statistical smearing due to the finite number of measured 
particles per event (|138p . Figure lllb compares the measured event-by-event 
distributions of V2 , v 3 and V4 measured by ATLAS (j!38|) , the distribution of the 
£ n from the IP-Glasma model (|7ip . and the final v n distributions after the viscous 
hydrodynamical evolution. The theoretical v n distributions after hydrodynamic 
evolution are seen to be slightly wider than the initial £ n distributions, indicating 
non-linear hydrodynamic effects on the fluctuation spectrum. These appear to 
be somewhat stronger in peripheral than in central collisions (|7ip . which may 
explain the stronger increase of a V2 / V2 compared to a £2 / £2 seen in Figure lllb 
for large impact parameters. In Figure lllb non-linear mode-mixing effects on 
the shape of the event-by-event distribution are particularly prominent for v^; 
only after accounting for the non-linear hydrodynamic evolution, the measured 
v n distributions can be fully reproduced. 

The ability of the IP-Glasma initial-state fluctuation model, which is deeply 
rooted in QCD, combined with viscous hydrodynamic evolution to simultane- 
ously correctly describe the mean ^2,3,4 coefficients at both top RHIC and LHC 
energies, their centrality and px dependence, as well as the detailed shape of 
their event-by-event distributions, is impressive. The recent measurements of 
these observables and their successful theoretical explanation has established the 
hydrodynamic paradigm for ultra-relativistic heavy-ion collisions beyond reason- 

3 Here we ignored flow angle fluctuations by setting v n = («„). 
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able doubt. Future measurements of the centrality and beam energy dependence 
of the complete v n power spectrum and of the correlations between the flow 
angles together with more complete theoretical modeling of the early-pre- 
equilibrium and late freeze-out stages, will allow for a precision determination of 
the transport properties of the QGP and the final completion of the Little Bang 
Standard Model. 

7 SYNOPSIS AND FUTURE PERSPECTIVE 

Construction work on the Little Bang Standard Model is nearing completion. 
Heavy-ion collisions generate many different classes of Little Bangs - the initial 
conditions depend on collision system, collision energy and collision centrality. 
Each class is characterized by its own spectrum of initial-state density fluctu- 
ations which (at LHC and top RHIC energy) appear to be calculable directly 
from QCD, using gluon saturation ideas. This so-called Glasma model also pro- 
vides a description of the early pre-equilibrium evolution of the Little Bang, 
although the absence of color field quantum fluctuations in the present imple- 
mentations, and the resulting lack of thermalization, are weaknesses that need 
to be fixed - corresponding work is going on while this is being written. Af- 
ter a very short time of around 0.2— 0.4fm/c (for top RHIC and LHC energies), 
the Glasma pre-equilibrium dynamics can be matched to viscous fluid dynam- 
ics of an approximately thermalized QGP, which then describes the main part 
of the fireball evolution (the next 5— lOfm/c, depending on collision energy and 
geometry) until hadronization. After hadronization, the macroscopic hydrody- 
namic evolution should be matched to a microscopic kinetic description based on 
hadronic degrees of freedom which then describes the rest of the Little Bang's life 
until final decoupling. At LHC energies, replacing this microscopic kinetic ap- 
proach by a continued application of fluid dynamics, with appropriately adjusted 
hadronic chemical potentials in the equation of state that describe the breaking 
of hadronic chemical equilibrium below T c , and with appropriately chosen values 
for the specific shear and bulk viscosities rj/s and £/s in the hadronic phase, is 
for many observables a reasonably good approximation, but will in the end not 
suffice for precision work. 

The fact that the initial Little Bang density profiles fluctuate from event to 
event is of crucial importance. Fluctuations in the magnitude and orientation 
of the initial-state eccentricity coefficients lead to fluctuations in the final har- 
monic flow coefficients and in the directions of those flows whose variances affect 
measured quantities in specific ways that require, for a faithful theoretical descrip- 
tion, event-by-event simulation of the hydrodynamic evolution. The correlations 
between the final flow planes differ from those between the initial eccentricity 
planes, sometimes dramatically, reflecting the non- linear hydrodynamic evolu- 
tion and mode-mixing, in particular with the large geometrically driven elliptic 
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flow in peripheral collisions. The recent measurement of these correlations and 
their successful explanation with full non-linear hydrodynamics, but not with 
linear response theory, constitutes an experimentum crucis in support of the hy- 
drodynamic paradigm. The community is looking forward to fully exploiting the 
spectrum of final anisotropic flow fluctuations and their angular correlations for 
pinning down the initial fluctuation spectrum and QGP transport properties. 

The present, still imperfect version of the Little Bang Standard Model is able 
to describe qualitatively, and in most cases even quantitatively, almost all ob- 
served features of soft hadron production in heavy-ion collisions at RHIC and 
LHC. The fact that collective flow anisotropics are uniquely sensitive the shear 
viscosity of the fireball matter has permitted the experimental determination of 
0?/ s )qgp ~ 0.2 = 2.5 x 7- at LHC temperatures (perhaps a bit smaller at RHIC 
temperatures), with less than 20% statistical and of order 50% systematic uncer- 
tainty related to the initial-state fluctuation and pre-equilibrium evolution mod- 
els. Rapid recent progress and ongoing work promise to reduce, within a year or 
two, the theoretical systematic error to a level competitive with the experimental 
statistical and systematic errors. Together with additional high-precision exper- 
iments, further exploration of the systematic variation of flow observables with 
collision energy, impact parameter and system size should permit the determina- 
tion of (?7/s)qgp with 5% relative precision, accurate enough to open the window 
for an experimental investigation of the (by an order of magnitude) smaller effects 
from bulk viscosity and the related shear and bulk viscous relaxation times. 
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Figure 1: Typical transverse energy density profiles e(x,y) from the IP-Glasma model (|8ip for a semipe- 
ripheral (& = 8fm) Au+Au collision at yfs = 200 A GeV, at times r = 0.01, 0.2, and 5.2fm/c. From 
r = 0.01 fm/c to 0.2 fm/c the fireball evolves out of equilibrium according to the Glasma model (175P76P771I78P : 
at r = 0.2 fm/c the energy momentum tensor from the IP-Glasma evolution is Landau-matched to ideal fluid 
form (for technical reasons (|7ip the viscous pressure components are set to zero at the matching time) and 
henceforth evolved with viscous Israel-Stewart fluid dynamics, assuming ?7/s = 0.12 for the specific shear 
viscosity. The pre-equilibrium Glasma evolution is seen to somewhat wash out the large initial energy den- 
sity fluctuations. The subsequent viscous hydrodynamic evolution further smoothes these fluctuations. The 
asymmetric pressure gradients due to the prominent dipole asymmetry in the initial state of this particular 
event (visible as a left-right asymmetry of the density profile in the left panel) is seen to generate a dipole 
("directed flow") component in the hydrodynamic flow pattern that pushes matter towards the right during 
the later evolution stages. 




Figure 2: Primordial fluctuation power spectrum of the Little Bangs created in 
2.76 A TeV Pb+Pb collisions of different centralities, from three different initial- 
state models (IP-Glasma, MC-Glauber, MC-KLN). 
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Figure 3: The eccentricity-scaled integrated elliptic flow of all charged hadrons, 
vf l (r)/s)/e, as a function of total charged hadron multiplicity density per unit 
overlap area, (1/ S){dN c \ 1 / dy). The experimental data points show two measures 
for the elliptic flow ((v 2 ) fl9Q]) and v 2 {2} (fTTOjl ) from 200AGeV Au-Au collisions 
at RHIC, measured by the STAR Collaboration. Both panels use the same sets of 
data, but use different normalization average initial eccentricities (e) and overlap 
areas {S) to normalize the vertical and horizontal axes, obtained from the Monte- 
Carlo-Glauber ([H]) and Monte-Carlo Kharzeev-Levin-Nardi (MC-KLN (]95tl93p ) 
models. The theoretical curves were computed with the VISHNU model (|17|,I1Q8[) . 
for different (temperature-independent) choices of the specific QGP shear viscos- 
ity (v/s)qgp- 
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Figure 4: Transverse mass (mr = ynrfi+pj) distributions of pions, kaons 
and protons from pp and Au+Au resp. Pb+Pb collisions at RHIC and LHC, 
at the collision energies indicated in the plots. The data (symbols) are from 
Refs. dHUIITlK (pp collisions at RHIC and LHC), Refs. (jngQIH) (Au+Au col- 
lisions at RHIC), and ([115]) (Pb+Pb collisions at the LHC). The models (lines) 
are described in the text. 
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Figure 5: Average transverse momentum vs. number of wounded nucleons N w , 
for pions, kaons and protons from 200j4GeV Au+Au collisions. Experimen- 
tal data from the PHENIX Collaboration are compared with (3+l)-dimensional 
event-by-event viscous hydrodynamic simulations with constant specific shear and 
bulk viscosities as indicated, for two different decoupling temperatures Tf = T^ in . 
Figure taken with permission from Ref. (|120p . 
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Figure 6: (a) Pion and (b) proton pt spectra from 200 A GeV Au+Au collisions 
at different collision centralities, compared with viscous hydrodynamic simula- 
tions (fTUSl) . 
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Figure 7: Two typical wounded nucleon configurations in the transverse plane 
for Au+Au collisions with N w = 100 wounded nucleons. Lines indicate contours 
of constant entropy density s after smearing the wounded nucleons with Gaussian 
distributions. Figure taken with permission from Ref. (I120p . 
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Figure 8: (a) Integrated elliptic flow at 2.76 TeV (|126|) in the 20-30% centrality class compared 
with results from lower energies taken at similar centralities. (b) The V2(pt) for pions and protons 
measured by STAR compared to hydrodynamic calculations with different eccentricities and rj/s 
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Figure 9: (a) The centrality dependence of v n {2} from 2.76 A TeV Pb+Pb collisions measured 
by ALICE (|102p compared to viscous hydrodynamic model calculations (|71|) . (b) Comparison of 
v u (pt) for the same collision system at 20—30% centrality from ATLAS (|133p with hydrodynamical 
calculations, using both a constant average and a temperature dependent rj/s (|7ip . 
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Figure 10: (a) The centrality dependence of the 5-particle cumulants from ALICE ()102p . 
The cumulant shown by green markers is sensitive to v^v^ cos[6( 1 I / 3— ^2)], the red markers to 
— v^vf cos[2(^2— ^1)]) and the blue markers to — V3V2V1 cos(3 l I / 3— 2^2 — ^l)- (b) Correlations be- 
tween different ^ n from ATLAS (|103p and from viscous hydrodynamics (jlOOp . Both data sets are 
from Pb+Pb collisions at the LHC. 
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Figure 11: (a) Relative event-by-event elliptic flow fluctuations versus collision centrality measured 
by ALICE (|145p . (b) Scaled distributions of t>2,3,4 (from top to bottom) from viscous hydrodynamics 
with IP-Glasma initial conditions (|7ip compared with experimental data from ATLAS (j 138|) and 
with the scaled distributions of the corresponding initial eccentricities £2,3,4- Nonlinear hydrodynamic 
evolution causes slightly larger variances for the v n distributions compared to those of e n . The data 
in both panels are from Pb+Pb collisions at the LHC. 



